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We design a system of phase oscillators that is able to produce temporally periodic sequences of patterns. 
Patterns are cluster partitions which encode information as phase differences between phase oscillators. The 
architecture of our system consists of a retrieval network with N globally coupled phase oscillators, a pacemaker 
that controls the sequence retrieval, and a set of patterns stored in the couplings between the pacemaker and 
the retrieval network. The system performs in analogy to a central pattern generator of neural networks and 
is very robust against perturbations in the retrieval process. 
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Sequential generation of patterns is an important 
issue in the context of dynamical systems that 
are designed for appHcations to artificial neural 
networks, gait models of animals, robot control- 
ling and the like. In this paper, we construct a 
system of phase oscillators that is able to gen- 
erate periodic sequences of patterns via imple- 
mented couplings and an inherent clock. The sys- 
tem consists of a pacemaker providing the clock, 
a retrieval network providing the patterns, and a 
set of couplings between the pacemaker and the 
retrieval network, coding the patterns. The pace- 
maker activates the stored patterns one by one as 
a function of its own phase, and the retrieval net- 
work selects the cluster partition corresponding 
to the activated pattern. The process of select- 
ing the cluster partition is realized by ensuring 
that the energy function of the retrieval network 
has only a single minimum while the pattern is 
activated. The gradient dynamics of this energy 
function then leads to a phase-locked motion of 
the retrieval network, such that the fixed phase 
relations between the oscillators represent the ac- 
tivated pattern. Thus the pattern amounts to a 
partition of oscillators into clusters characterized 
by a fixed phase relation (therefore called cluster 
partition). The system then performs similarly 
to a central pattern generator of biological neural 
networks. 



I. INTRODUCTION 

A sequential generation of patterns by a system of os- 
cillators is motivated by central pattern generators of bi- 
ological neural networks^i^. Central pattern generators 
usually are considered in the context of neuroscience as 
an explanation of how nervous systems produce move- 
ments. They are autonomous neural networks that can 



endogenously produce rhythmically patterned output like 
breathing, walking, or heartbeat. Apart from improv- 
ing an understanding of the biological aspects, nowadays 
the engineering aspect of networks of oscillators plays a 
prominent role^. For example, complex dynamical struc- 
tures in populations of phase oscillators were engineered 
by means of nonlinear time-delayed feedback that is im- 
plemented in the interactions of the oscillators'*. 

In this work we address the problem of generation of 
periodic sequences of patterns by a dynamical system of 
phase oscillators. We consider the case where a set of P 
patterns f should be presented periodically in a certain 
time order. To solve this problem we design a system 
which consists of a pacemaker, a retrieval network, and 
P stored patterns between these elements. The proposed 
system can be seen as an associative memory driven by 
a pacemaker. Here, we use a particular case of an asso- 
ciative memory (retrieval network) with time-dependent 
couplings controlled by the pacemaker through the pat- 
terns S,. Associative memory models based on the use of 
phase oscillators are alternatives to the Hopfield model^ 
with spins replaced by phase oscillators. These models 
are generalizations of the Kuramoto modcl^ (for a re- 
view see AcebronX). Examples for such generalizations 
were considered by Aoyagi^ and further improved by 
Nishikawa et al^ with the goal to achieve a storing ca- 
pacity similar to that of the Hopfield model along with 
error-free retrieval. 

As a result, our system is able to produce sequences of 
patterns in the retrieval network. These patterns are en- 
coded as phase differences between the oscillators. The 
associated dynamical states correspond to phase-locked 
synchronization. Such states amount to cluster partitions 
of the oscillators' set, in which oscillators sharing the 
same phase are gathered in one cluster. The system can 
be seen as a stylized version of the central pattern gener- 
ator of the leech heartbeat^, that is composed of two sets, 
the rhythm generator (corresponding to our pacemaker), 
and the pattern generator (the retrieval network). 

The paper is organized as follows. In section |TT] we 



present the model for our device and its different building 
blocks. In section IIIII we illustrate its performance and 
point on the analogy to a central pattern generator. The 
conclusions are drawn in section HVl 



II. THE MODEL OF OUR DEVICE 

The system is composed of a phase oscillator, which 
plays the role of pacemaker and is characterized by a 
phase variable ip, a retrieval network with N globally 
coupled phase oscillators with phases 4^, i = l,...,iV, 
and P patterns ^ stored in the couplings between the 
pacemaker and the retrieval network. 

The function of the pacemaker is to control the timing 
and activation of the stored patterns ^ as a function of 
its phase. In this way, the pacemaker works as a clock 
that points to the pattern that should be activated in 
the sequence. The retrieval network encodes the stored 
patterns ^ as phase differences between its N phase os- 
cillators, one by one as a function of time. We call it 
retrieval network since its function consists in retrieving 
the information stored in the couplings between the pace- 
maker and the network, and translating the patterns to 
the corresponding cluster partitions of the phase oscilla- 
tors. 

The operation of the retrieval network is achieved by 
designing a suitable energy function L. The minima of 
this function are approached by a gradient dynamics in 
which the energy function plays the role of the poten- 
tial. The minima correspond to specific cluster parti- 
tions, for which the phase differences take only values of 
zero or tt. These configurations are easily mapped to bi- 
nary sequences. In the set of N oscillators only A^ — 1 
phases are independent, therefore we can have 2^~^ clus- 
ter partitions. A pattern ^ is then described as an iV — 1- 
dimensional vector with components ^i being the phase 
differences between the oscillators i and N. 

Now, in order to retrieve a selected pattern out of the 
2^~^ ones that are in principle available, we implement 
couplings between the pacemaker and the phase oscil- 
lators which modulate the energy landscape in a way 
that the selected pattern becomes the only minimum at 
a given time. So the proposed system can be seen as an 
associative memory. The associative memory is particu- 
lar in the sense that the couplings are time-dependent, 
driven and controlled by the pacemaker through the 
choice of patterns ^. For the selection of another pat- 
tern, another set of couplings is addressed by the pace- 
maker. This way it becomes possible to retrieve a whole 
temporal sequence. 

If the retrieval network had only static sets of couplings 
between the oscillators, several minima could coexist and 
would be reached by the gradient dynamics depending 
on their basins of attraction. This case corresponds to 
an associative memory network as proposed by Aogagi^ 
and Nishikawa et al^. Our aim is, however, different. 
We consider time-dependent couplings in order to have 



a time-dependent energy landscape that is controlled by 
the pacemaker. In addition, we impose the constraint 
that only one minimum of L may exist at a time, when 
a pattern is selected by the pacemaker. The advantage 
then is that the gradient dynamics will lead to the re- 
quired minimum starting from any initial condition due 
to the absence of competing basins of attraction. 

In the following we shall describe the building blocks 
of our device in detail, before we illustrate how it works 
for retrieving a cyclic sequence of patterns. 



A. The retrieval network 

The dynamics of this network is characterized by an 
energy function L ("L" shall remind to its role as a Lya- 
punov function) whose gradient determines the phase 
evolution of the oscillators. In general, the system evolves 
to the different minima of L depending on the initial con- 
dition and on the basis of attraction of each minimum. 
For this retrieval network we use a particular version of 
the model proposed by Nishikawa et al^. We shall show 
how we tune the couplings in order to have only one min- 
imum in L when a pattern is selected for retrieval. The 
function L then depends on N oscillator phases $i in the 
following way: 
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with $j G [0, 27r[, i = 1, ..., iV the phase variables and K 
the coupling strength that finally determines the speed 
of convergence of the dynamics towards the stationary 
state, see Eq. [3] below. As noted before, out of the 
maximally N{N — l)/2 different phase differences, only 
A^ — 1 are independent, so we choose A$ijv = ^i — '^Ni 
i — 1, ...,N — 1 as independent variables, using the phase 
of the A^th oscillator as reference. With fij we denote 
the couplings, f^j e R, i,J £ {1, ..., A^}. 

For vanishing couplings fij we get back the Kuramoto 
model with twice the usual frequency in the interaction 
between the oscillators. Then, L has obviously 2^^^ min- 
ima given by the vectors S^'-''\k = 1,...,2^~^ of equal 
height with components Q e {0,7r}. In the configu- 
ration space of phase differences {A$iAr}, these minima 
are located at the corners of a hypercube of linear size tt 
with one corner in the origin of coordinates. It is these 
2^~^ local minima that are our candidates for retrieval. 
The label k of the state ^'^^^ is determined by the pattern 
of Os and tts interpreted as binary sequence in decimal 
representation. 

In general, the fij are couplings appropriately chosen 
to modulate the energy function L in order to retrieve 
a selected pattern ^^"^ A sufficient condition for a local 
minimum reads that the Hessian matrix of L is positive 
definite, i.e., all eigenvalues being larger than zero, due 
to an appropriate choice of fij . 



Now let us select one pattern ^('*) with s £ 
{1, ...,2^~^}. Let the couplings modulate the interac- 
tion between oscillator pairs (ij) according to 
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for i,j e {1, ..., N,i ^ j} with ^^' = 0, a any real num- 
ber with a > 1 and s the index of the selected pattern. 
So the couplings fij take values of ia, depending on the 
pattern. This choice is guided by the postulate that the 
selected configuration remains the only local minimum 
when these couplings are applied, while all other former 
local minima become saddles in at least one direction or 
local maxima. The conjecture is that a choice according 
to Eql2] satisfies this postulate. In the appendix we shall 
show that for this choice of couplings the Hessian is pos- 
itive definite for each configuration that is selected from 
the 2^~^ patterns, and not positive definite for all other 
2^^^ — 1 configurations, which were formerly also min- 
ima for vanishing couplings. (For a < 1 it can be shown 
that all local minima remain stable, but the selected one 
becomes the deepest.) The gradient dynamics will then 
retrieve the local minimum from any initial condition, 
not necessarily close to the selected minimum. 

It should be noticed that the choice of couplings fij for 
a selected pattern ^ by Eq. [2 can be mapped to the Heb- 
bian rule in the particular case where only one pattern 
is memorized by the associative memory of Nishikawa et 
al^. However, this restriction to the case of only one 
minimum of the energy function (ensured by the choice 
of a > 1) has the advantage that it leads to an error- 
free retrieval starting from an arbitrary initial condition, 
possibly far away from the final configuration. 

Explicitly the gradient dynamics of the oscillators then 
reads: 
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(3) 

Similarly to the Kuramoto dynamics the interaction of 
the oscillators depends only on phase differences A$y, 
and due to the choice of trigonometric functions the inter- 
action terms are bounded as in the Kuramoto dynamics. 
Due to the gradient dynamics the phase differences will 
evolve to a fixed point which is the minimum of the en- 
ergy function L that is closest to the initial conditions. 
Note that the second order Fourier term was related to 
the formation of clusters by MatoiS. In our case, this 
term appears due to the very construction of the energy 
function. 



B. The role of the pacemaker 

We extend our dynamics to sequential pattern retrieval 
via time dependent couplings controlled by a pacemaker. 



The pacemaker is a phase oscillator with constant fre- 
quency ajR and phase ip, whose time derivative is given 
by 



i> = wfl,. 



(4) 



It is the phase -ip that selects and activates a pattern 
r among the P stored ones over a duration B. This ac- 
tivation period lies in the phase (time) interval between 
Tp — Tpr and ip — ipr+i- For simplicity we have chosen 
iP^ = ^(r - 1) for ah r G {1, ..., P}, and B = ^. This 
means, when the instantaneous phase ip comes to the 
value V'r J couplings /^ are switched on that guarantee the 
retrieval of the pattern ''^'*(''). Here the pre-superscript 
r indicates the label of the pattern within the time se- 
quence, the post-superscript s{r) stands for the decimal 
label of the pattern r out of the selected subset. The 
corresponding dynamical equations read: 
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and 



r{^) = ei^p - ipr) - e(V' -i'r- B). 



(7) 



Here Q denotes the Heavyside function. This means 
that the function gr{i') controls the couplings to be given 

as a {^Y£,l — ""^^ I — 1) over the phase interval B, 

starting from ip = ipr on. 

Obviously we should ensure that ipr+i — "fpr ^ B > 
^Hrans, that is, the time interval between two initiations 
of pattern retrievals and the time of application of the 
constant couplings (to reach the new pattern) should be 
larger than the transient time iptrans which the retrieval 
network needs to go from one pattern to the next (the 
speed of the internal dynamics is controlled by K) . 

The information for the sequence generation is stored 
in the phase values ipr's through the functions 5r('0) 
which "switch on" the appropriate couplings. Since the 
stored patterns ^ are required at different times, they 
must be stored outside the retrieval network. There- 
fore we say that the patterns generated in the sequence 
are memorized in the couplings between the pacemaker 
and the retrieval network. This way the architecture can 
modulate the couplings fij between the oscillators of the 
retrieval network. 

In summary, the pacemaker transforms the static con- 
tents, memorized in the couplings, into a temporal se- 
quence of patterns, a feature that is in common with a 
central pattern generator. 

To complete our set of equations, we have to introduce 
noise in the retrieval dynamics. The reason is the fol- 
lowing. As we have seen, once a pattern is selected for 



retrieval, the system evolves according to the dynamics 
of Eq. [3] to the only stable fixed point of L of the re- 
trieval network. When the next pattern of the sequence 
is activated by the pacemaker, the system is still in the 
former fixed point that turns into an unstable one. To 
kick the system out of this fixed point and follow the re- 
quired sequence, we apply Gaussian white noise of small 
intensity T. Eq|3]is then replaced by 
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^^ = 'n ^ sin(A<D,,)(cos(A$^,)-/j,(^)j+Tr;,(t) 

(8) 
where rii{t) is a random variable describing the white 
noise with zero mean, < rii{t)rjj{t') >— 5ij5{t — t'), and 
T is the noise intensity. Our system is then described by 
Eq.smHH andH 



III. NUMERICAL STUDY 

In this section we study an example of this system 
and focus on its dynamical properties. We integrate the 
dynamics using a second order stochastic Runge-Kutta 
method^"'^ with time step Ai = 0.01 and noise intensity 
T = 0.001. We consider a set of A^ = 11 phase oscilla- 
tors in the retrieval network. This network can therefore 
encode 2^° = 1024 patterns as cluster partitions. Out 
of this set, we have randomly selected P — 5 patterns. 
They are ^e^^, 2^°, ''e^^, ^C^'^ '?'"'' in the indi- 
cated order. The time dependent-couplings are changed 
at phase values ipr — 0, 2tt/5, Att/5, 67r/5 and 8tt/5. 
The parameters are chosen a.s K = 10, B = 2tt/P, a ~ 2 
and ujR = 1. 

Figure [T^ shows the energy L as function of time eval- 
uated in the actual state of the retrieval network. L 
jumps from its minimal value at L w —500 to some larger 
value around L w —300, where it remains as long as the 
system of oscillators searches the new minimum, corre- 
sponding to the new choice of external fields. When the 
new minimum is found, L drops to the minimal value 
again. During such an interval of duration B, the Eu- 
clidean distance D{t) in configuration space between the 
actual state and the closest pattern (corner of the hyper- 
cube) has a peak at an intermediate time interval where 
the system is moving from one to the next selected state. 
We see these peaks in Fig. [TJs. For about half of the pe- 
riod B this distance is zero, indicating that the state of 
the system system corresponds to the required pattern. 

In Fig. [TJ; we plot the states which are closest to the 
instantaneous states of the system as a function of time. 
Obviously the closest states are just the selected ones, but 
this does not mean that the actual states (evolving with 
time) are identical with the selected ones over the whole 
duration of the plateau; as mentioned before, the distance 
to the selected states vanishes only for roughly half of 
the period as it is seen from Fig. [T|d. The width of the 
peaks in the distance from the closest states can be tuned 



by the coupling parameter K, large K accelerates the 
convergence to a new pattern, once the time-dependent 
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FIG. 1. Sequential pattern retrieval of a sequence si — {^^''^ , 
f, ?'*'*^ C'*^^ C^°^^}. (a) Energy L{t) as function of time. 

(b) Euclidean distance of the actual state of the system to 
the closest state (corresponding to patterns on the corners 
of the hypercube) . It vanishes for roughly half of the period 
B, so that the system then has retrieved the desired state. 

(c) Closest states to the current evolving state as function of 
time. As seen from the figure, the closest states themselves 
vary with time. The set of closest states agrees with the set 
of selected states. The states carry their decimal labels. 



In Fig. [2] we show three permutations of the same set of 
the five stored patterns. Indicated are the five plateaus 
in time where a certain pattern remains the closest to 
the current state and where this pattern agrees with the 
system's state over roughly half of the period (the analo- 
gous figures to Fig. [T}) are not displayed here) . It should 
be noticed that a change in the pattern sequence from 
Fig. [2k to [2}d and [21: only amounts to reorder the phase 
shifts ipr, no other change of the system's structure is 
needed. 
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FIG. 2. Sequential pattern retrieval of different sequences 
with the same states, (a) Time evolution of the sequence S2 
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(b) Time evolution of the se- 
(c) Time evolution of 



quence S3 

the sequence S4 = {^ 1023 ^^672^ ^942^ ^0^ ^477-^^ -^^ g^jj ^^^^ ^^ 

show only the states which are closest to the current evolving 
state. For roughly half a period the distance between the clos- 
est state and the actual state of the system vanishes, which 
is interpreted as pattern retrieval. 



IV. CONCLUSIONS 

We have designed a system of phase oscillators that is 
able to produce periodic sequences of patterns. Patterns 
are stored in the couplings of the system and retrieved 
and encoded as phase differences. Due to the task divi- 
sion between the pacemaker, the stored patterns and the 
retrieval network, the system is very flexible and robust. 
Different sequences of the stored patterns can be imple- 
mented without modifying the system's structure. The 
retrieval network itself operates in a robust way since it 
has by construction only one minimum, therefore the dy- 
namics converges to the desired pattern independently of 
the initial conditions. 

Our device may be regarded as a very stylized version 
of the central pattern generator of the leeches' heartbeat. 
According to Hooper—, the central pattern generator of 
the leech heartbeat can be divided into two sets, the 



rhythm generator (corresponding to our pacemaker) , and 
the pattern generator (corresponding to our retrieval net- 
work). The pattern generator there generates the actual 
motor pattern in response to the driving input from the 
rhythm generator (in our case in response to the driving 
input of time-dependent couplings from the pacemaker). 



V. ACKNOWLEDGMENTS 

We would like to thank Alexander S. Mikhailov very 
much for valuable discussions. 



VI. APPENDIX 

Consider the energy function for an all-to-all coupled 
system of N phase oscillators: 
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From now on we set the coupling strength K 
the couplings be chosen according to 
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^^ = xn = 0, a any real number 
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with a > 1 and s the index of the selected pattern. We 
now prove a sufficient condition that the Hessian matrix 
with respect to the N — 1 independent phase differences 
is positive definite for the selected pattern and not posi- 
tive definite for all other 2^~^ — 1 patterns, provided we 
choose the external fields /^ according to EqO From 
the first derivative dLjdxi we immediately see that can- 
didates for extrema are xi G {0, tt}, where Xi was defined 
as A$iAr = $i — $Ar, while for /i^r > 1 the individual 
cos-dependent terms are different from zero. For a par- 
ticular given choice of jij , fiN with possibly alternating 
signs the first derivatives can vanish also at intermediate 
values of Xi which we project on [0, 2tt[ that can lead to 
further extrema. This part we treat numerically in order 
to exclude that these extrema compete with the selected 
minimum that shall be retrieved. 

Next let us consider the Hessian of L as function of the 
phase differences. Apart from the normalization factor, 
its diagonal elements are given as: 

—J = -sin^(a:;fe) + cos(a;fc)( cos(a;fc) - fkN) 



dx 
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2J f sin^(a;fc - Xj) - cos{xk - Xj) (11) 
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its off-diagonal elements are: 



and in view of that we shall need a > 1, see below.) 



- — -— = ain^ixk - xi) - cos{xk - xi)(cos{xk - xi) - fik) ■ 
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(12) 
For a choice of the external fields /y according to the 
rule (HI, the matrix simplifies to 
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for the diagonal elements and fc = l,...,A^— 1, and to 
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for the off-diagonal elements, with k^l G {1,...,A^ — 
l},fc ^ I. The upper (lower) sign in front of the the 
bracket with fik stands for the case that the difference of 
components |a;| — xf\, read off from the selected configu- 
ration ^*, is zero (tt), respectively. 

Next we study the positive definiteness of this matrix 
for all 2^~^ configurations which may be selected for 
retrieval. Here it is convenient to classify the configu- 
rations in terms of their Hamming distance from the 
selected pattern, i.e. the number of mismatches of 
components between ^ and ^^ , which varies between zero 
and N — 1. By a suitable permutation of the coordinate 
axis in configuration space we can always achieve that 
the k mismatches occur in the first k coordinates of 
^ so that the corresponding Hessian H is chosen as 
representative for all patterns with k mismatches. 

H with no mismatches According to our choice 
of /ij, their signs are opposite to those of cos(A<I'^j,), 
that is A$fj. = or (tt), so that cosA$fj, = 1 or (-1) 
and fik = —a (or +a), a > 1, respectively. The diagonal 
elements then simplify to d^L/dx^ — {N — 1) (1 + a) , the 
off-diagonal elements to d^L/dxkdxi — — (1 + a). The 
Hessian therefore takes the form of an {N — 1) x (N — 1) 
dimensional circulant matrix, whose eigenvalues turn out 
to be Ai = 1 + a with multiplicity 1 and A2 = N{1 + a) 
with multiplicity (N-2). (Here we have used the follow- 
ing: Eigenvalues of an 71 x n circulant matrix, specified 
by the vector (cq, ci, ..., c„_i), are known to be given as 






rj27Tijk/n 



Ck with j = 0,-1,-2, 



^(n-1). 



In our case the Hessian has a particularly simple form, 
for which one element in each row is {N — 1)(1 + a), 
while all other N — 2 elements are —(1 -I- a). Using 
these values and the fact that the sum over all n roots 
of the unit circle adds up to zero leads to our results for 
the eigenvalues.) Now, since for a > 1 all eigenvalues 
are positive, the selected configuration corresponds to a 
local minimum in configuration space, whatever pattern 
has been chosen for retrieval. (In order to have only a 
local minimum at the selected configuration, obviously 
a > — 1 would be sufficient, but at the same time, the 
other patterns should become saddles or local maxima, 



H with one mismatch Next we evaluate the Hessian 
for a configuration that differs from the selected pattern 
in a single phase difference. Without loss of generality we 
assume the mismatch to happen in the first coordinate, 
affecting the Hessian in the first column and the first row 
according to Hu = {N —1){1 — a), Hij = (a — 1) = Hji 
for j = 2, ..., N -1, while the remaining {N - 2){N - 2) 
submatrix remains circulant. The Sylvester criterion, 
applied to the positive definiteness of the overall 
{N — 1) X {N — 1) matrix, is now violated due to the 
first element Hu = [N — 1)(1 — a) < for a > 1, so 
that the configuration with one mismatch is no longer a 
local minimum of the energy function L. (As necessary 
and sufficient condition for a Hermitian matrix to be 
positive definite, the Sylvester criterion requires that all 
leading principal minors of the matrix are positive.) 

H with fc > 1 mismatches Now the configura- 
tion has k mismatches with the selected configuration 
which we arrange to occur in the first k coordinates. 
Here it should be noticed that /^jv will have the "wrong" 
sign with respect to A$iAr, i = 1, ..,fc, but fu will have 
the "right" sign with respect to A^u for i, I e {1, ..., fc}, 
since two mismatches compensate in the relative phase 
differences ("wrong" (or "right") refer to the feature 
which prevents (or ensures) the property of becoming 
a local minimum, respectively.) This explains why the 
components of the fc x fc submatrix Sk (H) in the upper 
left corner of the Hessian are given by 

Sk{n) = {N~l)-a{N ~2k+l), i^l,...,k (15) 

for the diagonal elements and 

S'fe(Jj) = -(1 + a) i,j^l,...,k, ,iy^j (16) 

for the off-diagonal elements. The submatrix Sk{H) is 
again circulant and has eigenvalues Ai = {N — fc)(l — a) 
with multiplicity 1 and A2 = iV — a{N — 2k) with multi- 
plicity fc — 1, so that the determinant of this submatrix 
reads \Sk{H)\ — AiAj"^. Now we have to distinguish the 
following cases: 

1. k odd. For k odd, Aj"^ is always positive while Ai < 
for a > 1, so that \Sk{H)\ < for odd k and a > 1 and 
the Sylvester criterion for H being positive definite is 
violated as it should be for any positive number of mis- 
matches. 

2. k even. For k even, both eigenvalues may be negative 
so that \Sk{H)\ > 0. In order to see that the Sylvester 
criterion is still violated, we have to distinguish the fol- 
lowing cases: 

(i) For a > 1 and fc > N/2 we have Ai < and A2 > 0, 
so that the Sylvester criterion is violated, 
(ii) For a > 1 and fc < N/2, A2 < for a > N/{N - 2fc), 
so that \Sk{H)\ < only for 1 < a < N/{N - 2k). 
(iii) To finally see what happens for a > 1 and a > 
N/{N — 2k) let us consider the determinant of the sub- 
matrix of size Z = fc— 1 in the upper left corner of H. This 



matrix has eigenvalues cti = {N — fc + 1) — a{N — fc — 1) 
and (72 = N — a{N — 2k) with even algebraic multiplicity 
{k — 2), so that again the sign of Ai determines the sign 
of this subdeterminant. Now cti < for 1 < j^~^j]j^ < a, 
but this is certainly satisfied, since in the considered case 
k > 2 and a was even larger than N / [N — 2k) by assump- 
tion. So this (fc — 1) X (fc — l)-dimensional subdeterminant 
violates the Sylvester criterion for H to be positive defi- 
nite. 

In particular, for the maximal number of mismatches 
k = iV-1, Ai =: l-a < Ofora > 1 and A2 ^ N + a{N~- 
2) > ior N > 2, and for A^ = 2, X^'^ = A^"^ = 1 > 0, 
so that the corresponding pattern again ceases to be a 
local minimum of the energy function. 
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